Interplay between ferromagnetism and superconductivity in tunneling currents 



O 

o 

(N 
^■ 

00 



M. S. Gr0nsleth,i J. Linder/ J.-M. B0rven/ and A. Sudb0i 

'Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 

(Dated: Received February 2, 2008) 

We study tunneling currents in a model consisting of two non-unitary ferromagnetic spin-triplet 
superconductors separated by a thin insulating layer. We find a novel interplay between ferro- 
magnetism and superconductivity, manifested in the Josephson efi'ect. This offers the possibility 
of tuning dissipationless currents of charge and spin in a well-defined manner by adjusting the 
magnetization direction on either side of the junction. 

PACS numbers: 74.20.Rp, 74.50.+r, 74.20.-z 
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The coexistence of ferromagnetism (FM) and super- 
conductivity (SC) has recently been experimentally con- 
firmed 0, Qj 13- This offers the possibility of observ- 
ing new and interesting physical effects in transport of 
spin and charge. The spin-singlet character of Cooper 
pairs in conventional superconductors suggests that FM 
and SC are mutually excluding properties for a mate- 
rial. Indeed, coexistence of magnetic order with singlet 
superconductivity is not possible for uniform order pa- 
rameters 4] but can occur if the superconductor is in a 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) 5] state with 
modulated ferromagnetic and superconducting order pa- 
rameters characterized by a non-zero wave vector q. On 
the other hand, spin-triplet Cooper pairs i3| are in 
principle perfectly compatible with ferromagnetism. For 
instance, odd-in-frequency Sz = spin-triplet supercon- 
ductivity in superconductor-ferromagnet structures has 
been much studied in the literature 8]. The synthesis 
of superconductors exhibiting ferromagnetism, with si- 
multaneously broken U(l) and 0(3) symmetries, are of 
considerable interest from a fundamental physics point 
of view, and moreover opens up a vista to a plethora of 
novel applications. This has been the subject of theo- 
retical research in e.g. Refs. and has a broad range 
of possible applications. Also, focus on hybrid systems 
of ferromagnets and superconductors has arisen from the 
aspiration of utilizing the spin of the electron as a bi- 
nary variable in device applications. This has led to spin 
current induced magnetization switching |lO|. and sug- 
gestions have been made for devices such as spin-torque 
transistors and spin-batteries Moreover, spin 

supercurrents have a long tradition in '^He , while re- 
cent work has focused on dissipationless spin-currents in 
unitary spin-triplet superconductors |l4l |. 

This Letter addresses the case of two p-wave super- 
conductors arising out of a ferromagnetic metallic state, 
separated by a tunneling junction. Such states have 
been suggested to exist on experimental grounds, in com- 
pounds such as RuSr2GdCu2 08 UGe2 0, ZrZn2 
and URhCe Jj, and have been studied theoretically 
in e.g. Refs. [3 [13, E3- Coexisting FM and spin- 
triplet SC have also been proposed to arise out of half- 



metallic ferromagnetic materials such as Cr02, and the 
alloys UNiSn and NiMnSb . We compute the Joseph- 
son contribution to the tunneling currents, both in the 
charge- and spin-channel, within linear response theory 
using the Kubo formula. Our assumption is that the su- 
perconducting order is that of spin-triplet pairing, and 
we consider the analog of the A2-phase in '^He, i.e. SC 
order parameters that satisfy |Ak-f|| ^ |Akj,xl 7^ and 
Akti = 0. In terms of the dk- vector formalism _13i], 
we then have a non-unitary state since the average spin 
(Sk) = idk X d* = i(|AkTTp - |AkiJ')z of the Cooper 
pairs is nonzero. Such a scenario is compatible with 
uniform FM and SC since the electrons responsible for 
ferromagnetism below the Curie temperature Tm con- 
dense into Cooper pairs with magnetic moments aligned 
with the magnetization below the critical temperature 
Tc. The choice of such a non- unitary state is motivated 
by the fact that there is strong reason to believe that 
the correct pairing symmetries in the ferromagnetic su- 
perconductors (FMSC) discovered so far are non-unitary 
|l7l E(l Ell |. The exchange field will also give rise to 
a Zeeman-splitting between the t,! conduction bands, 
thus suppressing the SC order parameter Akj^ (3) as 
illustrated in Fig. QJ)). 

Another important issue to address is whether the SC 
and FM order parameters coexist uniformly, or if they are 
phase-separated. One possibility is that a spontaneously 
formed vortex lattice due to the internal magnetization 
m is realized in a spin-triplet FMSC [i^. However, there 
have also been reports of uniform superconducting phases 
in spin-triplet FMSC A key variable determining 

whether a vortex lattice appears or not is the strength of 
the internal magnetization m '24|. Current experimen- 
tal data on ZrZn2 and URhGe apparently do not settle 
this issue unambiguously, while uniform coexistence of 
FM and SC appear to have been experimentally verified 
in UGe2 113 • Furthermore, a bulk Meissner state in the 
FMSC RuSr2GdCu208 has been reported in Ref. [13, 
indicating the existence of uniform FM and SC as a bulk 
effect. Consequently, we will use bulk values for the order 
parameters and assume that they coexist uniformly. We 
emphasize that one in general should take into account 



2 



the possible suppression of the SC order parameter in 
the vicinity of the tunneUng interface due to the forma- 
tion of midgap surface states which occur for certain 
orientations of the SC gap. The pair-breaking effect of 
these states in unconventional superconductors has been 
studied in e.g. Ref. A sizeable formation of such 

states would suppress the Josephson current, although it 
is nonvanishing in the general case. Also, we use bulk 
uniform magnetic order parameters, as in Ref. j^^. The 
latter is justified on the grounds that a ferromagnet with 
a planar order parameter is mathematically isomorphic 
to an s-wave superconductor, where the use of bulk val- 
ues for the order parameter right up to the interface is 
a good approximation due to the lack of midgap surface 
states. Moreover, we consider thin film FMSC such that 
the Lorentz-force acting on the electrons will be unable to 
accelerate particles in a direction parallel to the junction. 
Our model is illustrated in Fig. ^). 

The main result of this Letter is that the Joseph- 
son current in spin- and charge-sector between two non- 
unitary FMSC can be controlled by adjusting the rela- 
tive magnetization orientation on each side of the junc- 
tion provided that spin-triplet Cooper pairs are present. 
Our system consists of two FMSC separated by an in- 
sulating layer such that the total Hamiltonian can be 
written as H = Hi^ + Hji + Ht, where L and R rep- 
resents the individual FMSC on each side of the tun- 
neling junction, and describes tunneling of particles 
through the insulating layer separating the two pieces of 
bulk material. The FMSC Hamiltonian is given by [l^ 
ffpMSC = Ho + EkV'k-^kV'k, where Ho = JN-f{0)m^ + 
^LkCT^kiT + LkQ/s ^kQ/3^kQ/3- Here, k is the electron 
momentum, Ekg. = Ek — o'Cz, o" =T:J.= ilj J is a spin 
coupling constant, 7(0) is the number of nearest lat- 
tice neighbors, m = {mxTrny^mz] is the magnetiza- 
tion vector, while Aka/3 is the superconducting order 
parameter and foka^a = (c-k^Cka) is the two-particle 
operator expectation value. The ferromagnetic order 
parameter is defined by C = 2J7(0)(rn2; — vmy) and 
C2 = 2J7(0)m2. We express the Hamiltonian in the basis 

V'k = (ckT Cki cLk| cLkj^)^, where Cko- (c^g.) are annihila- 
tion (creation) fermion operators. Note that there is no 
spin-orbit coupling in our model, i.e. inversion symmetry 
is not broken. Consider now the 4x4 matrix 

. _ 1 /-£kl + cr-C idk -cro-y \ 
2Uidk-Ta,)t (ea-^-C)V 

which is valid for a FMSC with arbitrary magnetization. 
As explained in the introduction, we will study in detail 
a non-unitary equal-spin pairing (ESP) FMSC as illus- 
trated in Fig. nji), i.e. Ak^x = Akxt = 0, C = in Eq. 

Since the quantization axes of the two FMSC are 
not aligned, one needs to include the Wigner d-function 
[30| denoted by 2?cr'o-(^?) with j = 1/2 to account for 
the fact that a f spin on one side of the junction is 



not the same as a t spin on the other side of the junc- 
tion. The spin quantization axes are taken along the 
direction of the magnetization on each side, so that the 
angle d is defined by mir • rriL = mpmL cos(t?) where 
rui — \m.i\. The tunneling Hamiltonian then reads ifx — 

Ekp..' 'D.'amT^^clj^a, +rkp4..CkJ, where we ne- 
gleet the possibility of spin-flips in the tunneling pro- 
cess. The validity of the tunneling Hamiltonian approach 
requires that the applied voltage across the junction is 
small. Here, we will be concerned with the case of zero 
bias voltage, so that the tunneling Hamiltonian approach 
is appropriate. Note that we distinguish between fermion 
operators on each side of the junction corresponding to 
Ckcr and dkcr- Furthermore, we write the superconducting 
order parameters as Akcrcr — lAkaale'^^''"'"^'"'-', where R 
(L) denotes the bulk superconducting phase on the right 
(left) side of the junction while 0k is a general (complex) 
internal phase factor originating from the specific form 
of the gap in k-space that ensures odd symmetry under 
inversion of momentum, i.e. 6^ = 0-k + tt. 




a) b) 



FIG. 1: a) Tunneling of Cooper pairs between two non-unitary 
FMSC with non-coUinear magnetization, b) Band-splitting 
for I, J, electrons in the presence of a magnetization in z- 
direction, leading to a suppression of interband-pairing. 

For our system, the Hamiltonian takes the form 
^^FMSC = Ho + Ha, Ha ^ X^k^ '/'L^k^'/'k^T, where we 
have chosen a convenient basis (j)^^^^ = (Ck^jC-ko-) that 
block-diagonalizes .Ak, and defined Ako- = \['^£\i.a<^z + 

Akcrcr (cTa; + iCj,) + ^\^„„{<7x — If J/)] with Pauli matriccs 
(Ti. This Hamiltonian is diagonalized by a 2 x 2 spin 
generalized unitary matrix C/km so that the supercon- 
ducting sector is expressed in the diagonal basis ^j^^ = 

4>i„U\^a = (7L'7_k<T)- Thus, Ha =^Eka <^L^ka^k,T, 
in which Akcr = U\^c,A-kaU^^ = diag(£;kcr, -£^ka)/2, and 
£^k. = V^L + IAkaaP. 

In order to find the spin and charge currents 
over the junction, consider first the generalized num- 
ber operator Nap = Sk'-ka'-k/S- The trans- 
port operator in the interaction picture then reads 

Ncpit) = -iEkp.[^-M'?)7kpcL(i)V(*)e""''' - 
^^a(i?)Tk*pdt^(i)Ck0(i)e""^], where eV = ^ll^- is an 
externally applied potential. The general current across 
the junction can be written 

l(t) = Y.Tap{Nap{t)), T = (-el, a) (2) 

a/3 
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such that the charge-current is l'^{t) = loit) while the 
spin-current reads = {Ii{t), l2it), l3{t)). Note that 
Eq. ||2Jl contains both the single-particle (sp) and two- 
particle (tp) contribution. The concept of a spin-current 
in this context refers to the rate at which the spin- vector 
S on one side of the junction changes as a result of tun- 
neling across the junction, i.e. S = i[iJT,S]. As there 
is no spin-orbit coupling in our system, this definition 
of the spin-current serves well [IJl ■ The spatial compo- 
nents of are defined with respect to the correspond- 
ing quantization axis. We compute the currents by the 
Kubo formula, (^a^(t)) = -i J^^^dt' {[N^fiit), HT{t')]), 
where the right hand side is the statistical expectation 
value in the unperturbed quantum state, i.e. when the 
two subsystems are not coupled. We now focus on the 
two-particle charge-current and z-componcnt of the spin- 
current, such that only a — (3 contributes in Eq. |(5J). 

Using linear response theory with the Mat- 
subara formalism, one arrives at {Naa{t))tp — 
2E^Sm{*a<T(eF)e-2'*'=^}, where ^^aAeV) is 
obtained by performing analytical continuation 
iLOn eV + iO+ {ujn = 2Tmk^T,n — 1,2,3...) on 
$„,(it^„) = -/„'/''^^^dre-"-(T,M„,(r)M,,(0)), where 
= EwpV^c.WTkpcijt)dp,it), while T is the 
temperature, and T,- denotes the time-ordering operator. 
Explicitly, we find 

kp 

X,p=± 

(3) 

where E'k^ = + l^kaaP, S,ka = £kff - Mr (R ^ L 

for k ^ p), and A^^^^(eF) = limi^„^ey+,o+ A[/(£^ka) - 
fiXpEp,)]/[iuj, + pE^a. - \Ep„]- X,p = ±1. 

In Eq. lO, we have used that T_k,-p = T^p which 
follows from time reversal symmetry, and f{x) is the 
Fermi distribution. Note that the chemical potential has 
been included in the excitation energies i^ka- In gen- 
eral, Eq. ^ will give rise to a term proportional to 
cos(0^^ —0^^), the quasiparticle interference term, in ad- 
dition to sin(6'^g. — 9^^), identified as the Josephson cur- 
rent. In the following, we shall focus on the latter while 
a coniprehensive treatment of all terms will be given in 
Ref. |33| ■ Consider now the case of zero externally ap- 
plied voltage {eV — 0). From Eq. we see that the 
Josephson charge-current becomes 

=e I] [1 + <Jacosd] cos((?p - ^k) sm{e^^ - O^J 

kpca 

X \Tkp\^\A^ac.\\^paa\F^p/{E^aEp^), (4) 

with F^; = E±[/(±^^ka) - f{Ep,)]/iEkc, T Ep^) while 
the expression for if ^{t) is equal except for a factor 
(— a/e) inside the summation. Observing that Eq. 
may be cast into the form /j^ = /q + Im cos(79), we have 
thus found a Josephson current, for both spin and charge. 



that can be tuned in a well-defined manner by adjusting 
the relative orientation i9 of the magnetization vectors. 
Below, we discuss the detection of such an effect. 

In the limit where one of the superconducting order 
parameters vanishes internally on both sides, i.e. the 
equivalent of an Al-phase, we see that the interplay be- 
tween ■& and ^CTCT remains, as only one term contributes 
to the spin sum over {ct, a}. In this case, the charge- and 
spin-current goes as cos^(i?/2) sin A^o-o-, where A9crcr = 
9^^ - 9^^ and Aq<^<^ with a = {T, i}, q = {k, p} is the 
surviving order parameter. For coUinear magnetization 
(i9 = 0), an ordinary Josephson effect driven by the su- 
perconducting phase occurs. Interestingly, one is able to 
tune this current to zero for niL || — niR (z9 = tt). 

Another result that can be extracted from Eq. Q is a 
persistent spin- Josephson current even if the magnetiza- 
tions on each side of the junction are of equal magnitude 
and coUinear — 0). This is quite different from the 
Josephson-like spin-current recently considered in ferro- 
magnetic metal junctions [29l |. There, a twist in the mag- 
netization across the junction is required to drive the 
spin-Josephson effect. In this Letter, however, we have 
found a persistent spin-current in the two-particle chan- 
nel even for coUinear magnetization. 

In the special case of eV = and equal SC phases on 
each side of the junction, i.e. 9^^ — 9^„, Eq. Q re- 
duces to the form /|_^ = Jq sin^(i?/2) sin(6'^ - 9^'.^) while 
= 0. This means that a two-particle spin-current 
without any charge- current can arise for non-coUinear 
magnetizations on each side of the junction in the ab- 
sence of an externally applied voltage and with equal SC 
phases 9]^^ = 9^„; see, however, Ref. 

It is well-known that for tunneling currents flowing 
in the presence of a magnetic field that is perpendicu- 
lar to the tunneling direction, the resulting flux thread- 
ing the junction leads to a Fraunhofer-like variation in 
the DC Josephson effect, given by a multiplicative fac- 
tor sin(7r$/$o)/(''''J*/^o) in the critical current. Here, 
$0 = T^h/e is the elementary flux quantum, and $ is the 
total flux threading the junction due to a magnetic field. 
However, this is not an issue in the present case since 
the magnetization is assumed to be oriented according 
to Fig. QJl). Since the motion of the Cooper-pairs is also 
restricted by the thin-film structure, there is no orbital 
effect from such a magnetization. 

Note that the interplay between ferromagnetism and 
superconductivity is manifest in the charge- as well as 
spin-currents, the former being readily measurable. Since 
the critical Josephson currents depend on the relative 
magnetization orientation, one is able to tune these cur- 
rents in a well-defined manner by varying d. This can 
be done by applying an external magnetic field in the 
plane of the FMSC. In the presence of a rotating mag- 
netic moment on either side of the junction, the Joseph- 
son currents will thus vary according to Eq. I^J. De- 
pending on the relative magnitudes of Iq and /„, the 
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sign of the critical current may change. Note that such 
a variation of the magnetization vectors must take place 
in an adiabatic manner so that the systems can be con- 
sidered to be in, or near, equilibrium at all times. Our 
predictions can thus be verified by measuring the critical 
current a.t eV — for different angles i9 and compare the 
results with our theory. Recently, it has been reported 
that a spin-triplet supercurrent, induced by Josephson 
tunneling between two s-wave superconductors across a 
ferromagnetic metallic contact, can be controlled by varv- 
ing the magnetization of the ferromagnetic contact [33j. 
Moreover, detection of induced spin-currents are chal- 
lenging, although recent studies suggest feasible methods 
of measuring such quantities |34l | . Observation of macro- 
scopic spin-currents in superconductors may also be pos- 
sible via angle resolved photo-emission experiments with 
circularly polarized photons 35| , or in spin-resolved neu- 
tron scattering experiments 36j. 

We briefly mention our results in the single-particle 
channel where we find that the charge-current and the 
z-component of the spin-current both vanish for eV = 0; 
see Ref. for details. They are nonzero for eV ^ 
even if the magnetization vectors are collinear. We 
stress that the finding of a non-persistent z-component 
of the spin-current does not confiict with the results of 
Rcf. j22j, as their z-direction corresponds to a vector 
in the xy-plane in our system. For Akcrcr 0, Isp(i) — 
ZEkpEa/j.^^.aW^^./sWlT^kppSml^^^A^-JC-eF)}, 
and the component of the spin-current parallel to 
riiL X niR is seen to vanish for i} = {0, tt} at eV = in 
agreement with Ref. [2^ . 

We reemphasize that the above ideas should be ex- 
perimentally realizable by e.g. utilizing various geome- 
tries in order to vary the demagnetization fields. One 
may also use exchange biasing to an anti-ferromagnet to 
achieve non-collinearity |37j . We have found an inter- 
play between FM and SC in the Josephson channel for 
charge- and spin-currents when considering non-unitary 
spin-triplet ESP FMSC with coexisting and uniform fer- 
romagnetic and superconducting order. 
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